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The exact 2-spinon part of the dynamic spin structure factor S xx (Q,ui) for the one-dimensional 
s=l/2 XXZ model at T=0 in the antiferromagnetically ordered phase is calculated using recent 
advances by Jimbo and Miwa in the algebraic analysis based on (infinite-dimensional) quantum 
group symmetries of this model and the related vertex models. The 2-spinon excitations form a 
2-parameter continuum consisting of two partly overlapping sheets in (Q, w)-space. The spectral 
threshold has a smooth maximum at the Brillouin zone boundary (Q = n/2) and a smooth minimum 
with a gap at the zone center (Q = 0) . The 2-spinon density of states has square-root divergences 
at the lower and upper continuum boundaries. For the 2-spinon transition rates, the two regimes 
< Q < Qk, (near the zone center) and Q K < Q < n/2 (near the zone boundary) must be 
distinguished, where Qk — > in the Heisenberg limit and Q K — » tt/2 in the Ising limit. In the regime 
Qk < Q < 7r/2, the 2-spinon transition rates relevant for S xx (Q,uj) are finite at the lower boundary 
of each sheet, decrease monotonically with increasing u, and approach zero linearly at the upper 
boundary. The resulting 2-spinon part of S xx (Q,ui) is then square-root divergent at the spectral 
threshold and vanishes in a square-root cusp at the upper boundary. In the regime < Q K < 7r/2, 
by contrast, the 2-spinon transition rates have a smooth maximum inside the continuum and vanish 
linearly at either boundary. In the associated 2-spinon line shapes of S xx (Q,uj), the linear cusps 
at the continuum boundaries are replaced by square-root cusps. Existing perturbation studies have 
been unable to capture the physics of the regime Q K < Q < n/2. However, their line shape 
predictions for the regime < Q < Q K are in good agreement with the new exact results if the 
anisotropy is very strong. For weak anisotropies, the exact line shapes are more asymmetric. 



I. INTRODUCTION 

Among all the spin-chain models that are directly rel- 
evant for the description of real quasi-one-dimensional 
(ID) magnetic insulators, the s = 1/2 XXZ model, 
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is the one whose physical properties have been studied 
most comprehensively. Today there exist more exact re- 
sults for this model than for any other model of compa- 
rable importance. n ™ 

The early demonstrationEJtl that the XXZ model is 
amenable to the Bethe ansatz led to a steady stream 
of advances in owp, understanding of many of its swoimd- 
state properties,Qi3 its thermodynamic pr|OpertiesJ3~El and 
the structure of its excitation spectrumBllj The T = 
phase diagram of the XXZ model, which was rigorously 
established by these advances, consists of a ferromagnetic 
phase at A > 1, a critical phase (spin- fluid, Luttinger 
liquid) at — 1 < A < 1, and an antiferromagnetic phase 
at A < -1. 

The mapping between the XXZ model and the exactly 
solvable 6-vertex and 8-vertex models yielded additional 



ground-state properties of the former on a rigorous basis, 
notably the spontaneous staggered magnetization in the 
antiferromagnetic phaae_and some critical exponents in 
the spin- fluid phase £3aO 

Until recently, exact results for the T — spin dynam- 
ics of the XXZ model were limited to a single nontrivial 
case, the XX model (A = 0). For this case, the spin sy^ 
tern is equivalent to a system of free lattice fermions,E3 
and the dynamic spin correlation functions—can be ex- 
pressed as fermion density correlati ons . (zz)Ej or as infi- 
nite determinants or Pfafnans (a;a:).r n 21 l In the surround- 
ing spin-fluid phase (—1 < A < 1), exact results for the 
infrared singularities of dynamic structiice factors were 
obtained by field-theoretic approaches Bt3 

A new avenue for the study of the T = dynamics 
of the XXZ model on a rigorous basis was opened up 
by important advances in the study of this model and 
the related vertex models in the framework of the alge- 
braic analysis based on quantum group symmetries. A 
detailed description of this method with all the results 
which our calculations build-on can be found in a recent 
book by Jimbo and Miwa.E3 Unlike the Bethe ansatz, 
this approach considers an infinite chain from the outset 
and exploits the higher symmetry of the infinite system 
(compared to the finite system) described by the quan- 



turn group U q (sl2)i 

The algebraic analysis of the XXZ model for the pur- 
pose of calculating correlation functions and transition 
rates (form factors) of local spin operators requires the 
execution of the following program: (i) span the infinite- 
dimensional physically relevant Hilbert space in the form 
of a separable Fock space of multiple spinon excitations 
and to generate the XX Z eigenvectors in this Fock space 
by products of spinon creation operators (so-called ver- 
tex operators) from the XX Z ground state (physical vac- 
uum); (ii) determine the spectral properties (energy, mo- 
mentum) of the spinon excitations; (iii) express the local 
spin operators in terms of vertex operators; and (iv) eval- 
uate matrix elements of products of vertex operators in 
this spinon eigenbasis. 

There exist two similar yet distinct programs which 
operate under different circumstances for essentially the 
same purpose. One is the fermion representation of 
the ID— s— j= 1/2 XY model or the equivalent 2D Ising 
modeljEiffij and the other is conforraal field theory for 
critical (massless) continuum models. Lj Quantum inverse 
scattering theory provides yet different ways of calcu- 
lating some correlation functions and matrix elements 
for massive r&lativistic continuum modelsO and for the 
XXX model.y 

The algebraic analysis^ operates in the massive phase 
stabilized by Neel long-range order at A < — 1, but 
the isotropic limit A — > — 1~ can be performed mean- 
ingfully at various stages of the calculation and thus 
yields equivalent-ijesults for the (massless) Heisenberg 
antiferromagnet .EJO 

In this paper we infer from the diverse ingredients now 
accessible via the Bethe ansatz and the algebraic analysis, 
an explicit expression for the exact 2-spinon part of the 
dynamic spin structure factor 



S xx (Q,lu) 



i oo 



(1.2) 



at T = and A < — 1. The line shapes thus ob- 
tained are of direct relevance for the interpretation of 
existing snectxpscopic data obtained via-inelastic neutron 
scatteringEao and Raman scattering^ on the quasi-lD 
magnetic compounds CSC0CI3 and CsCoBr3. 

In Sec. |lj we discuss the m-spinon eigenbasis and in- 
fer a suitable parametric representation of the energy- 
momentum relation for spinon excitations from it. In 
Sec. QI a closed-form expression for the 2-spinon den- 
sity of states D(Q,uj) is derived from this spectral in- 
formation. In Sec. IV we analyze the matrix elements 
(form factors) between the twofold degenerate ground 
state and the 2-spinon excitations and derive from them 
(in Sec. |vi|), after having solved the 2-spinon energy- 
momentum relations in the appropriate parametrization 
(Sec. M), a function M(Q,ui) which, when multiplied 
with D[Q, to), yields the 2-spinon dynamic structure fac- 
totS x 2 J(Q,u). 



A related study was previously undertaken by Weston 
and BougourziEl In that study the goal was to calculate 
the 2-spinon part of the dynamic structure factor S(Q, u)) 
defined as the Fourier transform of (cr n (t) ■ (To). The 
result was expressed as an expansion about the Ising limit 
(A — > —00) carried out explicitly to 12 th order. It is 
much more difficult to calculate this quantity than to 
calculate S X x{Q,oj), where no expansion is necessary to 
obtain explicit results. 

Finally, it is interesting to note that the exact result 
for the frequency-dependent spin autocorrelation func- 
tion <& xx (lu) = J^ 7T (dQ/2Tr)S xx (Q,uj) of the case A = 0. 
which was calculated in the fermion representation!!!! 
represents all m-spinon contributions for m — 2,4,... 
simultaneously. There the m-spinon structure of the ex- 
citation spectrum is reflected in $ I3; (w) by an infinite 
sequence of singularities at the band-edge frequencies 
u/J = 0,1,2,.... 



II. SPECTRUM 

The 2 Ar -dimensional Hilbert space of the XX Z model 
for a chain of N sites becomes non-separable in the limit 
N — ► 00. However, for the infinite chain, a separable 
subspace T can be constructed, and all physical prop- 
erties of the XXZ model can, in principle, be derived 
exactly from it. The classification of the XXZ spectrum 
in terms of m-spinon excitations, which is instrumental 
in the quantum group analysis, had_already been estab- 
lished by Faddeev and Takhtajan,E£lE3 for A = —1 in the 
framework of the algebraic Bethe ansatz. 

The (infinite-dimensional) space T is spanned by vec- 
tors |£ m ,e m ;.. . ;6.ei)j witn m = 0,1, ... and j = 0,1, 
which represent multiple spinon excitations. In the 
regime of interest here, the twofold degenerate vacuum 
state is represented by the two vectors |0)o, |0)i. These 
states break the translational symmetry of H . The trans- 
lation operator T (shift by one lattice site) transforms the 
two vectors into each other: 



T\0)j = |0) 



1-3) 



.7=0,1. 



(2.1) 



In the Ising limit (A — ► —00), they become the pure Neel 
states, |...UU ■■_■), I ...|TIT--->- 

Each spinon excitation is characterized by a (complex) 
spectral parameter £; and a spin orientation e; = ±1. The 
subspaces of T with even and odd numbers of spinon ex- 
citations are disconnected in all matters of concern here. 
They describe the physics of Apns with even and odd 
N asymptotically for TV — » oo.tiJ'EJ The. completeness re- 
lation for the spinon basis in T readso 
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J £i» ei)jj{£ij ei; ■ ■ ■ ; Cm, e m|- (2.2) 



These basis vectors are, in fact, eigenvectors of the XX Z 
Hamiltonian H and of the translation operator T 2 , 

m _. 

T\€ m ,e m ;...;t;i,ei)j = — — |£ m ,e m ;. . .;£i,ei)i_j, (2.3a) 



H\^ m £m,- ■ -;^l,^l)j — 2_^ e {£,i)\£m,£m- ■ ^iM) 



i = 1 



with the respective eigenvalues determined by 



r(0 = e 



-*p(0 



C 



-1 u 9" 



'(^" 2 ) ! 



.(O-ji^'logT®, 



(2.3b) 

(2.4a) 
(2.4b) 



in terms of the spectral parameter £ and the anisotropy 
parameter 



A = (q + q- l )/2, -Kq<0. 



(2.5) 



Here q is the.-deformation parameter of the quantum 
group U q (sl2)f-i and 



Qx (y) = (x;x)(y; x) (xy 1 ;x), 

oo 

(y;x)^l[(l-yx n ). 



(2.6a) 
(2.6b) 



n=0 



For most of the analysis to be carried out later, it is 
convenient to express £ in terms of the alternative spec- 
tral parameter /3: 



£ = ie 



in(3/2K 



2K </3 < 2K. 



(2.7) 



The energy and momentum of a spinon are then ex- 
pressed in terms of Jacobian elliptic functions, 



with 



e(0 = e(/3) = Jdn/J, 



r JK . , ttK' 
I = sinh . 

7T K 



(2.8a) 

(2.8b) 



(2.9) 



The anisotropy parameter (2.5) is related to the nomc 
-q = exp(-irK'/K), (2.10) 



and thus determines the moduli k,k' = Vl — k 2 of the 
elliptic integrals K = K(k),K' = K(k'). T he spinon 
energy-momentum relation resulting from (|2.8|), 



is equivalent 
Bethe ansatz 



eiO) =I\/l-k 2 cos 2 p, 0<p<n, (2.11) 

j corresponding relation obtained via 



For the calculation of SxJ (Q, to) from the 2-spinon den- 
sity of states and the 2-spinon matrix elements we intro- 
duce here translationally invariant vacuum states, 



10) 



|0)o + |0)i 
v/2 



|0)o-|0)i 
V2 



(2.12) 



which have wave numbers (total momenta mod 2n) 
and 7r, ,-sespectively, in the extended Brillouin zone 
(— 7r, +7t).Ej In the isotropic limit, the state |0) is a singlet 
(St = 0), and the state J7r) is the vector with S T = of a 
triplet (St = 1)- The corresponding linear combinations 
of 2-spinon states, 



i, /- n y _ |£2,e2;£i,ei)o + |£2,£2;£i,ei)i 
|?2,e2;?i,ei;U; = 



|£2,e2;£i,ei;i\ 



V2 
= |6 ; £2;Ci,ei)o" |6,e 2 ;£i,ei)i 

V2 



, (2.13) 



are then also translationally invariant, 



*b(«i)+p(6)] 



*b(?i)+p(«2 )+*■] 



|6,e2;£i I ei; 7r > 



r|£2,e 2 ;6,ei;0) = e 
r|£ 2 ,e 2 ;£i,ei;7r) = e 

Since the 2-spinon momenta and energies 

p(£i,6) = ^(/3i,/3 2 ) = Ma)+p(6), 

S(£i,6)-^(/3i,/3 2 ) = e(6) + e(6), 



|6,e 2 ;a,ei;0), (2.14) 



(2.15) 



are independent of the spin orientations ei,e 2 = ±1, all 
2-spinon states at fixed P will be at least fourfold degen- 
erate. In the isotropic limit, this degeneracy involves a 
singlet state (S T — St — 0) and the three vectors with 
S T = 0, ±1 of a triplet state (S T = 1). 

The four sets of 2-spinon excitations are readily identi- 
fied in the framework of the Bethe ansatz. In a finite sys- 
tem (N < oo), the singlet-triplet degeneracy is removed, 
and for anisotropic coupling (A < —1), the triplet levels 
are split up as well. The fourfold degeneracy emerges 
only asymptotically for N — ► oo, and thus reflects the 
higher U q (sl2) symmetry of the infinite system, which is 
used in the algebraic analysis. 



III. DENSITY OF STATES 

Here we c onsid er any one of the four sets of 2-spinon 
excitations (2.13) with fixed spin orientations £1,62 and 
express their energies ^(£1,^2) = s.i(pi) + e2(p2) in terms 
of the wave number Q = p\ + pi (0 < Q <j-37r) and the 
variable A = \(j>i - p 2 ) (-tt/2 < A < vr/2):EJ 



e 2 (Q, A) = ei(Q/2 - A) + ei(Q/2 + A). 



(3.1) 



These states form a continuum in (Q,c<j)-space, which is 
depicted in Fig. EL It consists of two partly overlapping 
sheets C± with boundaries 



wo(<2) 
w±(Q) 



2/ 

1 + K 

2/ 

1 + K 



■ sm ( 



yl + k 2 ± 2k cos Q, 



where 



k = cosQ K 



1-fc' 



(3.2a) 
(3.2b) 

(3.3) 



is the natural anisotropy parameter in most of the re- 
sults presented here. The excitation gap AE = 21k' 
approaches zero exponentially in the isotropic limit :□ 



AE ^ AnJexp 



(-tt 2 /v / -8(1 + A)) . (3.4) 



I 1 




FIG. 1. Two-spinon excitation spectrum ( |3.l| ) for k = 0.99 
(A ~ —2.305). It consists of two partly overlapping sheets 
C- and C+. Sheet C+ lies between u>o(Q) and ui-{Q) in the 
range Q K < Q < tv — Q K and between w+(Q) and u)-(Q) for 
vr — Q K < Q < n. Sheet C_ is obtained from C+ by reflection 
about the line Q = n/2. The interval (0, 7r) represents one 
half of the extended Brillouin zone, in which one ground state 
vector, |0), is assigned wave number Q = 0, and the other 
ground-state vector, \tt), wave number Q = n. 



The 2-spinon density of statescil 

tt/2 



D(Q,w) = ± J d\5(oj-e 2 (Q,X)) 



-7T/2 



£1 



was evaluated before in closed form:! 

D(Q,oj) = D+(Q,uj) + D_(Q,w), 
D±(Q,uj) = n±(Q,u>)/d±(Q,u>), 

for (Q,oj) 6 C±, respectively, and where 



(3.5) 



(3.6a) 
(3.6b) 



n±(Q,w) 

d±(Q,uj) 



2u 2 -(1 + k 2 )lu 2 ±2TcosQ 
4sin 2 Q ' 

T [(1 + k 2 )uj 2 -lu 2 {1 + cos 2 Q) T 2TcosQ] 1/2 
2 sin 2 Q ' 



T(Q,u) = Juj 2 - k 2 oj 2 Jll> 2 
With the auxiliary quantity 



W±(Q,uj) = J^k 2 - (?-±cosQ 



the result ( |3.6b ) can be written more compactly: 

w[sin 2 Q~WI(Q,uj)} 



D±(Q,u>) 



2T(Q,u>)W±(Q,w) 



(3.7) 



(3.8) 



(3.9) 



Note the reflection symmetry: D±(Q,lu) — D t (it—Q,uj). 
The 2-spinon density of state has square-root divergences 
all along the lower and upper boun daries of each sheet. 
At the zone center, expression (3.S) turns into 



D_(0,w) 



VIP 



Wu 



AI 2 k 12 



(3.10) 



The function D + (Q, u>) is plotted in Fig. for two values 
of anisotropy. 



IV. MATRIX ELEMENTS 

The m-spinon eigenbasis provides a useful framework 
for the separate analysis of the m-spinon contributions 
(to = 0,2,4,...) to any zero-temperature dynamical 
quantity of interest if a means of calculating the rele- 
vant matrix elements can be found. Here we focus on 
the 2-spinon matrix elements of the dynamic spin struc- 
ture fact or S XX {Q, m) at T = 0. 

With (|E3) the 2-spinon part of fll.2| ) has the form 

i +oo 

Si 2 l(Q,u)= 1 ~ E E E [dte^+W 

j,i'=0,ln=-ooei ) e 2 =± _ co 
2 

(4.1) 

The evaluation of this expression requires that we know 
all transition matrix elements of the spin operators o~ = 
i(ag ± zctq) between the vacuum states |0)j and the 2- 
spinon states \ti,2 , £2', £.1 , e i)j- All nonvanishing matrix el- 
ements of this type turn out to be related to each other, 



j<0|o3=|fc,T5fc,T>i = i-M°Z\&Mi,±)i-h (4.2a) 
j{$i,±\b,±Wt\Q)j=M<>t\-lZuT\-q&T)i, (4.2b) 




CO// 



FIG. 2. Normalized 2-spinon density of states D+(Q,ui) as 
a function of frequency for wave numbers Q K < Q < n and 
anisotropy parameter k — 0.7435 (A ~ —10) and k — 0.99 
(A~ -2.305). 



and can be expressed by a single function 
which was determined by Jimbo and Miwa:l 



(4.3) 



xHb,^) = & (I3 2 ,f3i) = ff 









where 



7(0 



tf=(«W) 

oo 

(x;y;z)= ][ (1 



(q^q 4 ^ 4 )^- 1 ^ 4 ^ 4 )' 
2(q 4 ;q 4 ;q 4 ) 



(<z 6 ;9 4 ;<z 4 )' 



xy n z m ) 



(4.4) 

(4.5a) 
(4.5b) 
(4.5c) 



Carrying out the space-time Fourier tran sform and the 
sum over the spin orientations in Eq. (4.1) yields 



i 2 r 2K r 2K 

x{d[Q + P(0i,0 2 )]\X o (^,^)+X\f32,0i)\ 2 

+5[Q-7r + P(/3 1 ,/? 2 )]|X°(/3 2 ,/3 1 )-X 1 ( / 3 2 , / 3 1 )| 2 }, (4.6) 



<->xa: (Q)< 



where we have also substituted ( |2.7| ). 



V. ENERGY-MOMENTUM RELATIONS 

Performing the integrals over the spectral parameters 
in expression (4.6) brings the 2-spinon dynamic structure 
factor into the form: 



4Sc«.»>-5EE3£S 



C=±<T=± 



J c (Q,u>) 



where the numerator, 



(5.1) 



(5.2) 



is governed by the 2-spinon transition rates and the de- 
nominator, 



'IK 

J C (Q,^) = 2( — 

1 7T 



dE dP 

W1W2 



dE dP 

W2W1 



pm 



21 



(2kK\ 






— 



|sn/3Jcn/3Jdn/3£ - sn/^cn/^dn/^l , (5.3) 



by the 2-spinon density of states. In these expressions, 
the spectral parameters now have fixed values {Px,l3$). 
These values are the solutions of the 2-spinon energy- 
momentum rel atio ns arising from the two products of 
S- functions in (fhq): 



u = E(p u 2 ), -Q = P(p 1 ,p 2 ) 
w = E{fc,fa), 7r-Q = P(^i,/3 2 ), 



(5.4a) 
(5.4b) 



for a = =F, respectively. Equations (5.4) with (2.15) and 
fl2.8|) are combined into 

u/I = dn/3i + dn/3 2 , (5.5a) 

— a sin Q = sn/?icn/3 2 + cn.j3isn.02, (5.5b) 

for the future analysis. 

For fixed a and at a generic point (Q,ui) within the 
range of the 2-spinon continuum, there exists exactly one 
distinct solution per sheet C±. Every such solution has 
multiplicity 8, acco unte d for by the permutation sym- 
metry f3\ «-> fa of ( |5.5| ) [factor 2] and the periodicity 
of the elliptic functions [fa ctor 4] . Now we use addition 
thcorcmscj to convert (5.5) into 



n,m— 



u 

1 



2dn/? + dn/3_ 

1 — /c 2 sn 2 /3+sn 2 /3_ 

. 2sn/3 + cn/3 + dn/3_ 

-asing= - - k2sn2 p +sa2 p_ 

with /3± = (/3i ± /3 2 )/2. From the ratio 
ct sin Q sn/? + cn/3 + 



w/J 



dn/3 4 




(5.6a) 
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(5.6b) 


3 2 




k i 


(5.7) 
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i-Q 
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(5.8a) 
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(5.8b) 


3 , 



These solutions yield 

dn/3+ = 



T 



T 



(5.9) 



and reduce (5.6a), effectively, into a quadratic equation 
for dn/3_ with c-independent solutions 



dn/3_ = 



1 ± cos Q u; 2 — klo 2 + T 



sinQ 



KLUn 



T' 



(5.10) 



where ujq(Q) and T(Q,uj) are given in ( p.2a| ) and (p.7[), 
respectiv ely. Finall y, a L anden transformation (k — * k) 
converts ( |5.9| ) and (p3n) into more explicit solutions in 
terms of incomplete elliptic integrals, 



0%{Q,u) 



1 + K 

-a — - — b 



. w 
arcsm — , k 

LU 



0°_(Q, U ) = ±±HF 



arcsm 



2IuW c 

k(1 + k)u>q 



(5.11a) 
, (5.11b) 



, and the new label c = ± 



where W C (Q, to) is given in (|3 
indicates that (Q,u>) 6 C±. 

The function (3~^(Q,uj) — (3Z(n — Q, u>\, which alone 
enters the final result, is plotted in Fig. || for two val- 
ues of anisotropy. It is finite along the lower bound- 
ary of C+ , decreases monotonically with increasing lu at 
fixed Qi an( A vanishes in a square-root cusp at the upper 
boundary. Note the different behavior along the por- 
tions u>+(Q) and u>q(Q) of the lower boundary of C_|_, 
P-(Q,u + ) = K,P^(Q,(jJ ) < K, which will give rise to 

different singularities in Sxx (Q, u>) in the two parts of the 
spectral threshold. 



VI. DYNAMIC STRUCTURE FACTOR 



k=0.7435 



O) 



CO. 





QA 



-a. 




(Si/I 



FIG. 3. Explicit solution [3^(Q,lu) of the 2-spinon en- 
ergy-momentum equations as functions of frequency for wave 
numbers Q K < Q < n and anisotropy parameters k = 0.7435 
(A ~ -10) and k = 0.99 (A ~ -2.305). 



as a product of a density-of-state function D±(Q, u>) and 
a transition-rate function M±(Q,uj), in generalization of 
the representation used in Ref. ^l| for the XXX model, 
and in analogy to the representation used in Ref. O for 
S zz {Q,uj) of the XY modehEl 



Si 2 x \Q^) = \Y, D c (Q,u)M e (Q,u). 



(6.1) 



With the solutions ( 5.11 ) of the energy-moment um rela- 
tions, the numerator (5.2) and the denominator (5.3) can 
be evaluated in the form: 



B c CT (Q,u/) = 



2K(k) 



-T]6 a 



l+ccosQ|A_(/3i)| 5 



A. Exact result for SiJ(Q,ui) 

J c (Q,u>) 
The 2-spinon dynamic structure factor ( fS.lj) for both 
sheets C± of the 2-spinon spectrum will now be evaluated respectively, where 



1— K. 

1 + K 



AK{k) 



-T]5 a _ 



UjTW r 



(6.2) 



(6.3) 






(6.4) 



with 7 = nK'/K, and i}d(% ) is a N evil le theta f unct ions 
With the exact results (3.9), (6.2), and (6.3), the 

physically motivated factorization (6.1) of SxJ(Q,u>) — 



5Z C =± S c (Q,ui) can now be established: 



D e (Q,u) = 



4K(k) 



M c (Q,u>)= ]T 

cr=±l 



2 , ,2 



-I 2 



4JC(k) 



sin 2 Q-W 2 {Q,lj)} 
2lu 2 J c (Q,oj) 



2lo 2 B°(Q,u;) 
; 2 [sm 2 Q-W 2 (Q,uj)}' 



(6.5) 



(6.6) 



The e xact 2-spinon transitio n ra te function thus obtained 
from (6^) substituted into (6.6) is 



M c {Q,lo) 



l - K 2 u$ + T 1 + ccosQ $ 2 A (pt) 
(1 + k)lj 2 sin 2 Q-W?# 2 d (P-)' ( ' ' 



where $ 2 A (/3) = \A-{[5)\ 2 with Q ) is the function 
& A {fi) =exp 



y^ e 11 cosh(2 7 ^) cosftTp - 1 , (( . s) 



- 1 I sinh(2/7) cosh(7?) 



(31 is given in (|5.11t|) , and t = 2/3/if. The final result 
for the exact 2-spinon dynamic structure factor reads: 



s?HQ,«)-^t 



? 9 2 



^(/31)|tan(Q/2)|- 



E 



(6.9) 



B. Line shapes and singularity structure 

The function M+(Q,u), which represents the 2-spinon 
transition rates for (Q,u>) 6 C + , is plotted in Fig. || for 
two values of anisotropy. The product of M + (Q, to) with 
the 2-spinon density of states D + {Q, to) (already shown in 
Fig. for the same two cases) yields the spectral- weight 
distribution S+(Q,w) of the 2-spinon dynamic structure 
factor for (Q,u>) <G C+. This function is plotted in Fig. ^ 
for four values of anisotropy, including the values chosen 
in Figs. §§ 

The transition rate function M + (Q,u>) exhibits qual- 
itatively different properties in the two regimes Q K < 
Q < n — Q K and n — Q K < Q < ir, where the spectral 
threshold is given by lj (Q) and lu + (Q), respectively (see 
Fig. fy). In the first regime, M + (Q,uj) is nonzero at the 
lower boundary, decreases monotonically with increasing 
w, and approaches zero linearly at the upper boundary. 
In the second regime, by contrast, M + (Q,to) approaches 




CO// 



FIG. 4. 2-spinon transition rates M+(Q, co) as a function of 
frequency for wave numbers Q K < Q < 7r and anisotropy pa- 
rameters k = 0.7435 (A ~ -10) and k = 0.99 (A ~ -2.305). 



zero linearly at both boundaries and has a smooth max- 
imum in between. The transition rates for (Q,u>) € C_ 
are the exact mirror image: M_(7r — Q, lj) = M + (Q,ui). 
These properties of M±(Q,lu) imply that the 2-spinon 

dynamic structure factor Sxx {Q, u>) diverges all along the 
portion u>o(Q) of the spectral threshold and that the lead- 
ing singularity there is the square-root divergence of the 
2-spinon density of states. Here the function S+(Q,lj) 
decreases monotonically from infinity at the lower bound- 
ary to zero at the upper boundary. Along the portion 
uj + (Q) of the lower boundary and along the entire up- 
per boundary of C+, the linear behavior of the transition 
rates removes the square-root divergence of the density 
of states in the product and replaces it by a square-root 
cusp in the dynamic structure factor. Here the spectral- 
weight distribution at fixed Q has a smooth maximum 
between the band edges. For strong anisotropy, the line 
shapes are broad and featureless. At moderate to weak 
anisotropy, the line shapes are distinctly asymmetric with 
the maximum positioned close to the spectral threshold. 
The function S x J(Q,lu), which is symmetric about 
Q = it/2, is equal to one or the other of the two func- 
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FIG. 5. 2-spinon dynamic structure factor SxJ (Q, oj) for (Q, u>) £ C+ as a function of frequency for wave numbers Q K < Q < ir 
and anisotropy parameters k = 0.7435 (A ~ -10), k = 0.8279 (A ~ -7), k = 0.8976 (A ~ -5), and k = 0.99 (A ~ -2.305). 



tions S±(Q,uj) except for Q K < Q < n — Q K . Here the 
two sheets overlap, and the 2-spinon dynamic structure 
factor has two square-root cusp singularities, one at the 
upper boundary of each sheet. The two upper boundaries 
coincide only at Q = ir/2. 



C. Isotropic limit 

When we take the isotropic limit A — > — 1~ in the re- 
sults for the 2-spinon density of states, transition rates, 
and dynamic structure factor for the purpose of link- 
ing up with results of calculations that were performed 
for the Heisenberg model in the (non-degenerate) critical 
ground-state, we must heed the fact that the size of the 
Brillouin zone changes from (— n/2, +tt/2) to (— n, +n) 
as the Neel long-range order in the ground state of the 
infinite system vanishes. 

In practice, this means that we switch our perspec- 



tive from considering both sheets C± of 2-spinon exci- 
tations over the range (— ir/2, +ir/2) to considering only 
the sheet C+ over the extended range (— it, +ir). The 
term with c = — in all expressions that contain a sum 
^ c=± is then be omitted. These contributions are now 
accounted for in the term with c = + over the extended 
range of physically distinct wave numbers. 

In the isotropic limit, the boundaries of C+ turn into 
the familiar sine curves, 



^o(<9) — > u>l(Q) = -^JsinQ, 

U-{Q) -> u>u(Q) = 7rJsin— , 

and the 2-spinon density of states becomes 

1 



D+(Q,u) 



y/LJ*(Q) - " 2 



(6.10) 



(6.11) 



A major simplification occurs in the 2-spinon matrix 
elements, 



x l (b,b) 



*°(6,6), 



(6.12) 



whic h implies that all terms with a = — in ( p.l[ ) and 
(3.2) disappear in the isotropic limit.. Th e results pre- 
sented previously for the isotropic caseMEil can be recov- 
ered from the more general result pr esen ted here, if we 
replace the spectral parameter (3 in (2.7) by the scaled 
spectral parameter a = —f3ir/2K', and then take the 
l imit q — > — 1 — . In this limit, the auxiliary expressions 
(4.5) can be simplified considerably: 



4U/4 



7(0 



(g 4 ;g 4 )(l-g 4 ) 

( g W) 2 r( 3 M+(W2) - 

(? 4 ;9 4 )(i 



q 4 ) 3 / 4 A-(a) 



where 



A±(a) = exp 



T{\)A-{m/2) 



dx sinh [x{\ + ia/ir) 
x sinh 2x cosh x 



,+•>' 



(6.13a) 
(6.13b) 

(6.14) 



represents (3.4) with x = jl in the limit 7 — » 0. 
The magnitude of the resulting 2-spinon matrix element 
becomesEj 



|*°(&,6) 



n 
2o" 



A-(ai - a 2 ) 



sinh (f + ^) sinh (f + %) 



Expression (6.7) for M + (Q,w) reduces to 



M + {Q. : uj) = -e-^\ 



where 



T f°° , cosh 2a; cos xt — 1 

/(£) = / aa; — - 1 o — — — e", 

Jo 



.r 



sinh 2x cosh x 



(6.15) 



(6.16) 



(6.17) 



.,/ I ,_a lu>l(Q)-ul(Q) 



— = -(«!- a 2 ) = cosh 



^ 



«£(Q) 



In Ref. [3l]this expression was evaluated, and its-unplica- 
tions for the 2-spinon dynamic structure factored 

S^(Q,u)=M+(Q,u)D+(Q,u), A = -l (6.19) 

were discussed in considerable detail. 

D. Ising limit 

When we analyze the exact 2-spinon dynamic struc- 
ture factor for very strong anisotropy (A — » —00), a 
convenient expansion para mete r about the limiting Ising 
model is K as defined in (|3.3[) . Expressed in terms of 



this parameter, the exchange anisotr opy ( |2.5| ) and the 
amplitude of the 2-spinon dispersion (2J3) become 



IAI 



;i + o( K 2 )) 



1 



(1 + k + 0{k 2 )) . (6.20) 



Here we set J = 1. To lowest order in k, the overlap 
region of the two sheets C± that make up the 2-spinon 
continuum (see Fig. uh collapses to a single point at Q — 
7r/2. The continuum boundaries, as expressed in terms 
of the reduced frequency 



n = w/2i-i, 

are now described by the curves 

Q±(Q) ~* ±kcosQ. 



(6.21) 



(6.22) 



The regime Q K < Q < it — Q K of the 2-spinon continuum 
is thus squeezed to zero width at Q = 7r/2, where the 
2-spinon bandwidth is zero as well in lowest order. 

The expansion of the 2-spinon dynami c st ructure fac- 
tor can be carried out in the final result (6^9) by using 



W c (Q,w) -»fi;sin 2 QWl 



n 



kcosQ 



MP-) 



MP- 



-2, ,2 



cos 2 Q ' 



2 sin (2/3f 



AluoWr 



(1 + k)kuJq ' 
yielding the closed-form expression 



S<$(Q,uj) 



1 



2 cos 2 Q 



lcos 2 Q- [ - 



(6.23) 



(6.24) 



(6.25) 



(6.26) 



which is identical (in lowest order) to the result 



( 6 - 18 ) S™(Q, U ) 



y/Acos 2 Q~{cj-\A\) 2 
4 cos 2 Q 



N ! 1_ ^t( cosQ + w ~I a I) 



(6.27) 



obtained by Ishimura and Shibat3 from a first-order per- 
turbation calculation about the Ising limit. 

It is instructive to perform the expa nsio n at an ear- 
lier stage of the calculation, namely in (L£), which then 
becomes 



Sg>(Q,w) 



1 



dpi 



d(i 2 s : m 2 {2(3-)5(uj-E) 



[S(Q + P) + S(Q - 7T + P)] 



(6.28) 



with 



E(pi,fo) -> - [l + K - K(sin 2 (3x + sin 2 #,)] 

and where we have used 

^°(/3i,/3 2 )^0, 
|X 1 (/3i,/3 2 )| 2 ^4sin 2 (2/3_). 



(6.29a) 
(6.29b) 



(6.30a) 
(6.30b) 



The asymptotic 2-spinon energy-momentum relations 
(5.5) are 



u + 2 + 2/k = 2(sin 2 fa + sin 2 /3 2 ), 
-crsinQ = sin(/3! + /3 2 ), 



and the solutions 
sin(2/^) 



-ersinQ, 



sin(2/3i) = Wl- 



n 



kcosQ 



(6.31a) 
(6.31b) 



(6.32a) 
(6.32b) 



have multiplicity 8 for (Q,u>) within the boundaries 
( 3.22 ) of the asy mpto tic 2-spinon continuum. Performing 
the int egrals in (6.28) then yields the asymptotic version 
of Q: 



S xx (Q,t 



E 



sin 2 (2/31) 



^2|sin(2 / 3i)cos(2/3^)| 



(6.33) 



from which (6.26) is recovered upon evaluation. 



E. Experiments 

Two of the most intensively studied physical real- 
izations of the ID s = 1/2 Heisenberg-Ising antiferro- 
magnet are the quasi-lD magnetic compounds CSC0CI3 
and CsCoBr3. A very comprehensive set of spectro- 
scopic data, which probe diverse aspects of the low- 
temperature spin dynamics of these materials, is now 
available from several inelastic neutron scattering ex- 
periments p|SKformed over the,-Gpurse of 15 years-at the 
BrookhavepEp, Chalk RiverBcJ, Laue-LangevinO, and 
RutherforcO Laboratories. This impressive collection of 
data is complemented by a set of high-resolution spec- 
troscopic data for the long- wavelength spin dynamics ob- 
tained via Raman scattering.EJ 

The basis for the interpretation of all the experimental 
data that involve the frequency range now known under 
the name 2-spinon continuum has be en t he perturbation 
calculation about the Ising limit q£,(L1) carried out to 
first order by Ishimura and Shibaliij (see also Ref. 47), 
which yielded the explicit expression ( 6.27 ) for the T = 
dynamic structure factor S xx {Q,lo). Whereas this cal- 
culation does reproduce the 2-spinon continuum bound- 
aries correctly to first order in the expansion parameter, 
the reliability of its line-shape prediction - a broad peak 



with the maximum near the center of the band and steep 
drops near both boundaries - has remained very much in 
question. 

The fact is that the line shapes observed in all exper- 
iments turned out to be highly asymmetric with a high 
concentration of intensity near the spectral threshold and 
a tail extending to the upper continuum boundary. 

Various attempts have been made at reconciling the 
discrepancy between theory and experiment h* consid- 
ering a second order perturbation calculatioro for the 
pure XX Z model, and by considering the impact .of bi- 
axial anisotropyEljaext-nearest-neighbor coupling)^] in- 
terchain coupling^ and exchange mixing,cJ all within 
the framework of a first-order perturbation treatment. 

In the range — 10 < A < — 7 of anisotropies, which 
best describes CSC0CI3 and CsCoBr3 according to some 
indicators, the line-shape predictions obtained via per- 
turbation calculation are in fair agreement with the ex- 
act 2-spinon result for wave numbers near the zone cen- 
ter. However, for wave numbers near the zone boundary, 
the exact 2-spinon spectral-weight distribution is more 
asymmetric, consistent with the experimental data. For 
Q K < Q < 7r — Q K , S xx (Q,uj) even diverges at the spec- 
tral threshold and decreases monotonically toward the 
upper continuum boundary. The conspicuous line-shape 
asymmetry found in the experimental data near the zone 
center, which is not at all reproduced by the perturba- 
tion calculation, does also exist in S xx '(Q,uj), but only 
for weaker exchange anisotropy (see Fig. pf). 
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